The long journey from the giant-monopole resonance to the 
nuclear-matter incompressibility 



J. Piekarewic^] 

Department of Physics, Florida State University, Tallahassee, FL 32306 

(Dated: February 8, 2008) 

Abstract 

Differences in the density dependence of the symmetry energy predicted by nonrelativistic and 
relativistic models are suggested, at least in part, as the culprit for the discrepancy in the val- 
ues of the compression modulus of symmetric nuclear matter extracted from the energy of the 
giant monopole resonance in 208 Pb. "Best-fit" relativistic models, with stiffer symmetry ener- 
gies than Skyrme interactions, consistently predict higher compression moduli than nonrelativistic 
approaches. Relativistic models with compression moduli in the physically acceptable range of 
K = 200 — 300 MeV are used to compute the distribution of isoscalar monopole strength in 208 Pb. 
When the symmetry energy is artificially softened in one of these models, in an attempt to simulate 
the symmetry energy of Skyrme interactions, a lower value for the compression modulus is indeed 
obtained. It is concluded that the proposed measurement of the neutron skin in 208 Pb, aimed 
at constraining the density dependence of the symmetry energy and recently correlated to the 
structure of neutron stars, will also become instrumental in the determination of the compression 
modulus of nuclear matter. 
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The compression modulus of symmetric nuclear matter is a fundamental property of the 
equation of state. While some of the existent claims in the literature may be overstated — 
indeed, there is little evidence in support of a correlation between the compression modulus 
and the physics of neutron stars [|1| — the compression modulus impacts on a diverse set 
of phenomena ranging from nuclear structure to supernova explosions. In particular, the 
compression modulus controls the energetics around the nuclear- matter saturation point. 
This is because the first derivative of the energy-per-nucleon with respect to the density 
(i.e., the pressure) vanishes at saturation, so the dynamics of small density fluctuations 
around the equilibrium position becomes solely determined by the compression modulus. 

To date, most efforts devoted to the study of the compression modulus have relied on 
the excitation of the isoscalar giant-monopole resonance (GMR). While the first set of mea- 
surements of the GMR date back to the late seventies and early eighties 0, [|, a recently 
improved a-scattering experiment finds the position of the giant monopole resonance in 
208 Pb at £ GM r = 14.17±0.28 MeV @. While the experimental story on the GMR in 208 Pb 
seems to be coming to an end, the theoretical picture remains unclear. On the one hand 
nonrelativistic calculations that reproduce the distribution of isoscalar-monopole strength 
using Hartree-Fock plus random-phase approximation (RPA) approaches with state-of-the- 
art Skyrme || ||] and Gogny |7| interactions, predict a nuclear compression modulus in the 
range of 7^ = 210 — 220 MeV . On the other hand, relativistic models that succeed in repro- 
ducing a large body of observables, including the excitation energy of the GMR, predict a 
larger value for the nuclear incompressibility (K ~ 275 MeV) [§, |9] . It is the aim of this 
paper to elucidate the origin of this apparent discrepancy. It is proposed that this discrep- 
ancy, at least in part, is due to the density dependence of the symmetry energy; a poorly 
known quantity that affects physics ranging from the neutron radius of heavy nuclei to the 
structure of neutron stars ||10| . It should be noted that while knowledge of the symmetry 
energy is at present incomplete, the proposed measurement of the neutron radius of 208 Pb 
at the Jefferson Laboratory |TT| should provide stringent constraints on this fundamental 
component of the equation of state. 

In this paper we follow closely the philosophy of Blaizot and collaborators who advocate 
a purely microscopic approach for the extraction of the compression modulus of nuclear 
matter from the energy of the giant-monopole resonance |7|, [12| . While the merit of macro- 
scopic (semi-empirical) formulas for obtaining qualitative information on the compression 
modulus is unquestionable 
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the field has attained a level of maturity that demands 
stricter standards: it is now expected that microscopic models predict simultaneously the 
compression modulus of nuclear matter as well as the distribution of isoscalar monopole 
strength. Moreover, theoretical studies based solely on macroscopic approaches have been 
proven inadequate 
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The starting point for the calculations is an interacting Lagrangian density of the following 
form: 



Ant = ^ [<?s0 - (<?v^ + • b M + ~(1 + T 3 )A^ 



A 



This Lagrangian includes an isodoublet nucleon field (ip) interacting via the exchange of 
scalar ((f)) and vector (V 1 , b M , and A 11 ) fields. It also incorporates scalar-meson self- 
interactions (k, and A) that are instrumental in reducing the unreasonably large value of 
the compression modulus predicted in the original (linear) Walecka model [17], |18 . The 



Lagrangian density depends on five unknown coupling constants that may be determined 
from a fit to ground-state observables. Four of these constants (g s , g v , k, and A) are sensitive 
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to isoscalar observables so they are determined from a fit to symmetric nuclear matter. The 
four nuclear bulk properties selected for the fit are as follows: i) the saturation density, ii) 
the binding energy per nucleon at saturation, Hi) the nucleon effective mass at saturation, 
and iv) the compression modulus (see Table |). It is noteworthy, yet little known, that the 
above four coupling constants can be determined algebraically and uniquely from these four 
empirical quantities . It is also possible for the various meson masses to enter as undeter- 



mined parameters. However, here the standard procedure of fixing the masses of the u and 
p mesons at their physical value is adopted; that is, m v = 783 MeV and m p = 763 MeV. As 
infinite nuclear matter is only sensitive to the ratio g^/m 2 , the mass of the a-meson must be 
determined from finite-nuclei properties; the cx-meson mass has been adjusted to reproduce 
the experimental root-mean-square (rms) charge radius of 208 Pb (r c h = 5.50±0.01 fm.) 

The symmetry energy of nuclear matter is a poorly known quantity with an uncontrolled 
density dependence in nonrelativistic models (for a recent discussion of the symmetry energy 
in Skyrme models see Refs. [^UJ and 0]). In contrast, the symmetry energy displays a weak 



model dependence in relativistic approaches. It is given by the following simple form: 

k 2 a 2 P 

s ^ = ink + ^r-2—2> ( 2 ) 



where E* F = \/k^ + M* 2 . The symmetry energy, together with its density dependence, is con- 
strained in relativistic approaches because the only "free" parameter in Eq. (g) is the NNp 
coupling constant. As the effective nucleon mass M* has been fixed in symmetric nuclear 
matter (and spin-orbit phenomenology demands a value in the range of M*/M = 0.6 — 0.7) 
reproducing the empirical value of the symmetry energy at saturation (J~37 MeV) con- 
strains the NNp coupling constant to a relatively small range. Note that relativistically, 
the density dependence of the symmetry energy can also be modified through the inclu- 
sion of isoscalar- isovector couplings terms [10[, density-dependent coupling constants J22 



and isovector-scalar mesons [23]. For simplicity, however, none of these contributions will 
be considered here. In reality, the symmetry energy at saturation is not well constrained 
experimentally. Rather, it is an average of the symmetry energy near saturation density 
and the surface symmetry energy that is constrained by the binding energy of nuclei. Thus, 



a prescription first outlined in Ref. |T0| is adopted here: the value of the NNp coupling 
constant is adjusted, unless otherwise noted, so that the symmetry energy at fcp = 1.15 fm -1 
(i.e., p = 0.10 fm~ 3 ) be equal to 26 MeV (see Table |). 

The nuclear observables used as input for the determination of the model parameters are 
listed in Table |. In all cases the saturation density, binding-energy-per-nucleon, and rms 
charge radius in 208 Pb have been fixed at their empirical values. Thus, the only discrimi- 
nating factors among the three "families" are the effective nucleon mass and the symmetry 
energy. While best-fit relativistic models suggest values for the symmetry energy and its 
slope at saturation density satisfying J > 35 MeV and L > 100 MeV, respectively |13 



family C is defined with an artificially small value for J (and correspondingly for L) in a 
"poor- man's" attempt at simulating nonrelativistic Skyrme forces |2l| . That nonrelativistic 
Skyrme models have a softer symmetry energy is revealed by the behavior of one of the 
most sensitive probes of the density dependence of symmetry energy: the neutron skin of 
208 Pb. Indeed, the neutron skin of 208 Pb is predicted to be equal to R n —R p = 0.16 fm for the 
recent SkX parametrization and falls below 0.22 fm for all eighteen Skyrme parameter sets 



considered in Ref. [20]. In contrast, best-fit relativistic models consistently predict larger 



values. For example, the NL3 model of Ref. ||, the TM1 model of Sugahara and Toki [24 
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and the NLC model of Serot and Walecka p5| , predict R n — R p = 0.28,0.27, and 0.26 fm, 
respectively (see also Table [TT]). 

Within each family defined in Table Q, calculations of the isoscalar monopole response 
have been performed using a compression modulus in the physically acceptable range of K — 
200—300 MeV. To illustrate the similarities and differences between these three families, the 
equation of state for symmetric nuclear matter (left panel) and the symmetry energy (right 
panel) are displayed in Fig. |T] at K = 250 MeV. Clearly, the properties of symmetric nuclear 
matter at saturation density are identical in all three models. Further, having fixed the 
value of the effective nucleon mass in symmetric nuclear matter, the full density dependence 
of the symmetry energy is determined by one sole number: its value at fcp = 1.15 fm -1 . 

Results for the peak energy of the giant-monopole-resonance in 208 Pb as a function of the 
nuclear incompressibility are listed in Table [TJ and displayed in Fig |2|. All calculations were 
performed using the nonspectral, relativistic random-phase-approximation (RPA) approach 
of Ref. [p6fl . For each family, there is a clear correlation between the compression modulus 
and the energy of the GMR. Indeed, all of the results are well represented (in this limited 
range of K) by a linear relation with a "universal" slope: 

£gmr = £200 + 0.026(7^ - 200) , (3) 

where Equr, E 2 oo, and K are all given in MeV. The intercept is non- universal and given 
by: E 2 oo = 12.22 MeV, £ 2 oo = 12.71 MeV, and E 200 = 13.14 MeV, for families A, B, and C, 
respectively. 

A few comments are now in order. First, the value of the slope (0.026) is obviously 
small. This suggests that even without theoretical uncertainties, it would not be possi- 
ble to determine the compression modulus from the 208 Pb measurement alone to better 
than A£ , G mr/0.026 MeV (AE GMR is the experimental uncertainty). At present, the best 
determination of the peak position of the GMR is E GMR = 14.17±0.28 MeV ||, thereby 
resulting in an uncertainty in the compression modulus of about 20 MeV. Second, and more 
importantly, the journey from the GMR to the compression modulus is plagued by uncer- 
tainties unrelated to the physics of symmetric nuclear matter. To illustrate this point we 
invoke, although never use in any of the calculations, a semi-empirical formula based on a 
leptodermous expansion of the nuclear incompressibility: 

K(A, I) = K + Ksutf/A 1 / 3 + K sym I 2 + K Co uiZ 2 /A 4 / 3 + . . . , (4) 

where K snr {, K sym , and i^coui are empirical surface, symmetry, and Coulomb coefficients 
and I = (N — Z)/A is the neutron-proton asymmetry. The sizable contribution from the 
surface term to K(A, I) has been discussed recently by Patra, Vinas, Centelles, and Del 



Estal [27| in the context of a relativistic Thomas-Fermi theory so we limit ourselves to only 
a few comments. A surface dependence is modeled here through a change in the value of the 
effective nucleon mass (surface properties are also sensitive to the a-meson mass but this 
value has been chosen to reproduce the rms charge radius of 208 Pb). As shown in Table |, 
family A uses an effective nucleon mass of M*/M = 0.6 while family B uses M*/M = 
0.7; all other input observables are identical. A larger M* generates a slightly compressed 
single-particle spectrum and a correspondingly smaller spin-orbit splitting. Consequences 
of this change in M* result in a larger intercept, as displayed in Fig. ^. Thus, compression 
moduli of approximately i^ = 275 MeV (for family A) and K = 250 MeV (for family B) are 
required to reproduce the experimental energy of the GMR. Further, if one incorporates the 
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experimental error into this analysis, one concludes that "best-fit" relativistic mean-field 
models are consistent with a compression modulus in the range J\~ = 245 — 285 MeV. 

We now turn to the central idea behind this work, namely, how our incomplete knowledge 
of the symmetry energy impacts on the the extraction of the compression modulus. Let us 
then start by considering two identical models, but with vastly different symmetry energies, 
that predict a compression modulus of K = 250 MeV. Further, for simplicity we assume 
that these two models have identical surface and Coulomb properties so only the first and 
third term in Eq. (|j) are relevant to this discussion. Both models attempt to reproduce the 
"experimentally" accessible quantity: 

K 208 = lim K(A, 7=0.212) = K + Ji" sym (0.212) 2 + . . . , (5) 

A— *oo 

defined as the compressibility of infinite nuclear matter at a neutron-proton asymmetry 
identical to that of 208 Pb (see Table [H]). The first model, having a very stiff symmetry 
energy (that is, K sym large and negative) reduces K(A, I) from its 1 = value of 250 MeV 
all the way down to, let us say, 200 MeV at I = 0.212. Comparing this prediction to the 
assumed experimental value of K 2 os = 225 MeV, it is concluded that the compression modulus 
of symmetric nuclear matter must be increased to J\~~275 MeV. The second model predicts 
a very soft symmetry energy. So unrealistically soft, let us assume, that it generates no shift 
in going from 1 = to 7 = 0.212 (i.e., K sym = 0). In this case, the compression modulus must 
then be reduced to K = 225 MeV to reproduce the experimentally determined value. Thus 
two models, originally identical as far as symmetric nuclear matter is concerned, disagree 
in their final values of the compression modulus due to an incomplete knowledge of the 
symmetry energy. While the situation depicted in Fig. ||] might not be as extreme, it does 
follow the trends suggested by the above discussion. Indeed, family C, with the softest 
symmetry energy, generates the largest intercept and consequently predicts the smallest 
compression modulus of the three families. 

In summary, the impact of the poorly known density dependence of the symmetry energy 
on the extraction of the compression modulus of nuclear matter from the energy of the 
giant-monopole resonance in 208 Pb was addressed. The nuclear matter equation of state 
and the distribution of isoscalar monopole strength in 208 Pb were computed using three 
different families of relativistic models constrained to reproduce a variety of ground-state 
observables. For each family the compression modulus was allowed to vary within the 
physically acceptable range of 1^ = 200 — 300 MeV. The first family (A), with an effective 
nucleon mass fixed at M*/M = 0.6 is, at least for K = 275 MeV, practically indistinguishable 
from the successful NL3 model of Ref. M. The second family (B) differs from the first in 
that the effective nucleon mass is increased to M*/M = 0.7, thereby generating a slightly 
compressed single-particle spectrum but still a robust phenomenology. Finally, the third 
family (C) is obtained from the second one by artificially softening the symmetry energy in 
a "poor-man's" attempt at simulating nonrelativistic Skyrme models. When the peak energy 
of the GMR is plotted against the compression modulus, a linear relation with a universal 
slope is obtained. In contrast, the intercept is family dependent and it is largest for the model 
with the softest symmetry energy. Demanding agreement with the experimental value for the 
peak energy fixes the compression modulus at: i^ = 275, 255, and 240 MeV, for families A, B, 
and, C, respectively. It is therefore suggested that the discrepancy between relativistic and 
nonrelativistic models in the prediction of the compression modulus of nuclear matter may, 
at least in part, be due to our incomplete knowledge of the symmetry energy. At present, 
this issue can not be resolved. Yet the proposed Parity Radius Experiment (PREX) at the 
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Jefferson Laboratory should provide a unique constraint on the density dependence of the 
symmetry energy through a measurement of the neutron skin of 208 Pb. Such a measurement 
could have far-reaching implications: from the determination of a fundamental parameter 
of the equation of state (K) to the structure of neutron stars p0|. 
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TABLE I: Empirical bulk observables used in the determination of the coupling constants and 
the scalar mass. The symmetry energy J has been fixed at &f = 1.15 fm -1 but the quantities in 
parenthesis represent its value at saturation density. 



Family 


k° (fm" 1 ) 


e (MeV) 


M*/M 


K (MeV) 


J (MeV) 


r ch (fm) 


A 
B 

C 


1.30 
1.30 
1.30 


-16 
-16 
-16 


0.6 
0.7 
0.7 


200 
200 
200 


-300 
-300 
-300 




26(38) 
26(37) 
20(28) 


5.50±0.01 
5.50±0.01 
5.50±0.01 


TABLE II: The compression modulus of symmetric nuclear matter, the slope of the symmetry 
energy at saturation density, the compression modulus for asymmetric (I = 0.212) nuclear matter, 
the neutron skin of 208 Pb, and the energy of the GMR in 208 Pb for the three families discussed in 
the text. 


Family 


K (MeV) 


L (MeV) 


K 208 (MeV) 


Rn 


-Rp (fm) 


E GMR (MeV) 


A 


200 
225 
250 
275 
300 


120 
120 
119 
119 
119 


184 
203 
224 
246 
268 






0.28 
0.28 
0.28 
0.28 
0.28 




12.27 
12.88 
13.58 
14.14 
14.81 


B 


200 
225 
250 
275 
300 


108 
108 
108 
108 
107 


187 
208 
230 
252 
276 






0.25 
0.25 
0.26 
0.26 
0.26 




12.65 
13.35 
14.03 
14.75 
15.36 


C 


200 
225 
250 
275 
300 


82 
82 
82 
82 
82 


190 

212 
235 
258 
282 






0.19 
0.19 
0.19 
0.19 
0.19 




13.13 
13.80 
14.45 
15.09 
15.81 
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FIG. 1: Equation of state for symmetric nuclear matter (left panel) and the symmetry energy 
(right panel) as a function of the Fermi momentum for the three families discussed in the text. In 
all the cases presented here the compression modulus was fixed at K = 250 MeV. 
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FIG. 2: Energy of the isoscalar giant-monopole resonance as a function of the nuclear matter com- 
pression modulus for the three families discussed in the text. The box displays the experimentally 
allowed range of E GMR = 14.17±0.28 MeV @. 
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